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ABSTRACT 


<  /VI  ^ 


In  this  paper  the  properties  of  the  generalized  Euler-Frobenius 


polynomial 


are  studied.  It  is  proved  that  its  zeroes  are  separated 


- ^  ii  v  j  tr¬ 
oy  a  factor  q  and  their  asymptotic  behavior  as  (g  •»  °»j  is  obtained.  As  a  \ 

consequence  it  is  shown  that  least  squares  spline  approximation  on  a 

biinfinite  geometric  mesh  is  boundable  independently  of  the  (local)  mesh 

ratio  q  and  that  the  norm  of  the  inverse  of  the  corresponding  B-spline  Gram 


matrix  decreases  monotonly  to  2k  -  1  for  large  q,  as  /q+ 
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•  A  - 


SIGNIFICANCE  AND  EXPLANATION 


The  Euler-Frobenius  polynomial  II  plays  an  important  role  in  the 
analysis  of  cardinal  polynomial  splines.  It  has  simple,  negative  zeroes,  and 
the  fact 

W"1*  * 0 

allows  ua  to  conclude  that  there  exists  a  unique  bounded  cardinal  polynomial 
spline  f  e  S  that  interpolates  prescribed  data  b  e  l  at  integers. 

The  essential  properties  of  cardinal  polynomial  splines  have  been  later 
extended  to  the  more  general  case  of  cardinal  £-splines,  and  thereby,  by  an 
appropriate  change  of  variables,  to  polynomial  splines  on  the  biinfinite 
geometric  mesh 

(q1)^,  1  <  q  <  “  • 

The  generalized  Euler-Frobenius  polynomial  is  in  the  latter  case  given  by 


nn  ( X ;  q) 


.  n 
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nltn  i®0 
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n-if-n 


(")  I  I  (q 
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X ) ,  t  s  =  inq 


In  the  report  some  new  characteristics  of  H^( • ; q )  are  outlined.  A 
simple  but  far  reaching  property  is  proved:  the  zeroes  are  separated  by  a 
factor  q.  This  fact  helps  us  to  analyse  spline  interpolation  on  the 
biinfinite  geometric  mesh.  In  particular,  it  is  proved  that  the  least  squares 
spline  approximation  is  boundable  independently  of  the  local  mesh  ratio  q, 
and  that  the  norm  of  the  inverse  of  the  corresponding  B-spline  Gram  matrix 
decreases  monotonly  to  2k  -  1  for  large  q,  as  q  ♦  00 . 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive- - 

summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


ON  THE  GENERALIZED  EULER-FROBENIUS  POLYNOMIAL 
Y.  Y.  Feng*  and  J.  Kozak* 


1 .  Introduction 

The  exponential  Euler  polynomial  An(xit)  played  an  important  role  in  the  analysis  of 
cardinal  polynomial  splines.  This  is  much  due  to  the  fact  that  the  spline  defined  by  the 
functional  relation 

$  (x)  s-  A  (x»X>,  x  e  [0, 1  [  , 
n  n 

<p  (x  +  1)  tm  (x)  otherwise  , 
n  n 

vanishes  at  all  integers  for  particular  values  of  X,  the  zeroes  of  the  Euler-Frobenius 

polynomial  II  (X)  (1  -  X)nA  (OjX).  A  nice  survey  of  cardinal  polynomial  splines  can  be 

n  n 

found  in  [7].  Micchelli  (6}  showed  that  the  essential  properties  of  cardinal  polynomial 
splines  can  be  extended  to  the  more  general  case  of  cardinal  jC-splines.  By  applying  his 
results  to  the  particular  differential  operator. 


X  (D)  i-  TT  (D  -  it),  D  s-  —  ,  t  e  R 
t  l -0  dx 


and  to  the  corresponding  generalized  exponential  Euler  polynomial 


A  (xjX.q)  !-  -1—  l  (-)""1  (")  -3 -  ,  q  :■  e* 

"  nit"  i-0  1  q1  -  X 


he  analyzed  spline  interpolation  at  knots  or  the  biinfinite  geometric  mesh 

(g1)^  . 

In  this  case,  the  generalized  Euler-Frobenius  polynomial  is  given  by 
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2.  The  Zeroes  of  I!  (X;q) 
n 

We  start  the  section  with  the  symmetries  of  the  generalized  Euler-Frobenius 
polynomial.  In  addition  to  the  description  (1.3),  we  shall  use 


n-1  .  . 

I  an  ,(<5)^  - -  n  (X;q),  Y  :»  - - 

i-0  n,i  Y  (q  -  1)"  "  n  nit" 


to  emphasize  its  polynomial  character  in  X. 

Theorem  2.1.  The  polynomial  II  (Xjq)  satisfies 
"  n 


„  ,,  .  ,n-1  -n{n-1)/2_  ,  n, -1 

n  (Xjq)  -  X  q  n  (q  X  >q) 

n  n 


The  coefficients  an  ^ ( q )  can  be  recurrently  computed  by 

an+1,i<q>  “  (<l  -  “  <in”^)an>i<d)  +  <qn+W  -  Dan.i-^q))  .  (2.2) 


^  •n.O^1  *n,-1<q)  -  *n,n<q)  8“  0  * 

Proof.  For  n  «  1  or  X  «  0,  (2.1)  obviously  holds.  Assume  X  *  0,  n  >  2.  Then 


nn<X,q>  -  Y  I  <-)n“i(")  M  <qJ  -  X) 
i-0  j«0 


Yn  l  (->"'i(")q1X_1  I  I  (qJ  -  X)  , 


since  the  n-th  order  finite  difference  of  a  constant  vanishes.  But 


i.-l  ,  j  ,  ,  .  ,n.n-t  -n(n-1)/2  ,  j  n,-1, 

q  X  |  I  (qJ  -  X)  »  (-)  X  q  II  (qJ  -  q  X  ) 


which  completes  the  proof  of  (2.1). 


-.<«« i  i  jii'ii.i  wmf  H^pei 


T  -nr  ;  --yrr 


— .... 


Wliwfeifii 


In  terms  of  the  the  recurrence  relation  (1.4)  reads 

JQ  an+1,i<<1,xi  “  “(q_1  >_1((1-»qn  -  (q"+1-X)  a^^q)!*) 

-  (q-1)"1  I  ((qn+1-qn-1)a  (q)  +  (q"41-1-!).  (q))!1  , 

±mQ  n»1  nfl*l 


if  we  define  an>_,(q)  “  an<n(q)  8“  0,  and  thia  confirms  (2.2). 
Corollary  2.1.  The  coefficients  an  i(q)  satisfy 

Vi(«>  -  ‘JB<n'21"1)/2«„,„.i.1(q»  » 


and  for  n  >  2 


.n  l(q)  -  q<n-i><n-1-i>/2  a(J)qJ 


The  integer  coefficients  a  are  sytmnetric 

*  » . 


In  particular. 


C  -  ("I1)  * 


-  <„ .  2,(?)  - (£)  -  (JJ)  . 

It  is  easy  to  prove  (2.3)-(2.6)  by  using  (2.2)  and  mathematical  induction.  We  shall 
omit  this  step. 

Prom  now  on  we  think  of  the  seroes  of  ( • »q)  as  functions  of  q.  It  is  proved  in 
(61  that  the  n  -  1  seroes  Mn>i(q),  i  -  1,...,n  -  1,  of  ly-jq)  are  all  simple  and 
real,  in  fact  negative.  They  satisfy 

“n,!^’  «  °’  <  0  •  <2- 


lim  u  .(q)  -  0,  lim  u  fq)  -  all  i  . 


MWCTV;*  '  r«vr«"?' 


(2.9) 


y  (q  1 )  -  y  ^  . (q) ,  all  i  • 

n,l  n,n*-i 


We  ahall  think  of  the  y  . <g)  as  ordered, 

n.i 


.<<!>  K  U„  -><q)  <  ***  <  U„  „_1(q>  <  ° 

n,1  n,2  n,n-i 


Then,  additionally,  by  (4)  and  (2. IS) 


A  f  A  Wn  1(<}> 

4-  —  <  o,  4-  >  o 

dq  (  q  *1  _n_1 


(2.11) 


The  symmetry  (2.9)  tells  us  that  we  can  restrict  our  discussion  to  the  case  q  >  1. 
Lemma  2.1.  Let  q  >  1 .  Then 

yn  i_1(q>  <  wlvH,i,q)  <  qwn,i<q)'  1  “  2,3 . .  “  1»  B  ”  2'3 . 

Proof.  Suppose  y  .  , (q)  <  qu  ,(q)  holds  for  some  n.  By  hypothesis  then 

r—  n/i-i  Hpi 

sigmdl  (a  ^/q))  •  sign(II  (Xjq))  <0,  X  *  [qy  .(q),y  .  (q)l  , 

n  n  n  f  x  n  1 1 


(2.12) 


and  from  (1.4) 


nn+1(X,q)  *  °'  *  «  IqMnti(q),yn>1<q)] 


(2.13) 


Also  by  ( 1 .4) 


But  y  (q)  is  a  zero  of  II  (*;q),  thus  another  look  at  (1.4)  tells  us 
n  t  i  n 

sign(Hn+1  (qun  1(q>fq>)  *  Jn+1  (,1n  ,q> ,q)  >  <  0  ' 

and  there  is  at  least  one  zero  of  ( • iq)  in  each  of  the  intervals 

]y„  4  .(q) «qy  ,(q>[  •  ail  i  .  (2.14) 

n,i-i  n#i 

Also  by  ( 1 .4) 

signdl  ^.(O+iq))  •  sign (II  (y _ ,(q)»q)>  <  0  , 

n+i  n+i  Tifii"! 

sign  (IT  (y  (q);q)>  •  sign(H ,  (-*iq))  <  0  , 
n+i  n# l  nti 

and  this  reveals  the  position  of  the  smallest  and  the  largest  zero  of  ( • »q) .  However, 

H  ,(*;q)  is  a  polynomial  of  degree  <  n  *  1,  and  in  each  of  the  intervals  (2.14)  there 
n+ 1 

is  exactly  one  zero,  y  . .  ,(q). 

IV* T  pi 

Now  (2.15)  brings  the  induction  hypothesis  to  the  next  level  and  (2.12)  is  proved 
since  it  obviously  holds  for  n  ■  2.  ■ 

It  is  easy  to  deduce  the  following  interesting  properties  of  y  (q). 

n  p  i 


(2.15) 


Corollary  2.2.  The  zeroes  M  . (q)  of  H  (*;q)  have  the  following  properties 
■ _  n,i  n 


Mn,l(q)  ‘  Pn.n-i(q)  =  q  '  811  1 


In  particular 


P2k,k(q)  "  -*  ' 


(2.16) 


and  for  1  <  [■ 


n  -  1  i 


,  .  .  n-i  .  .  „  1 

W_  ,  <q>  <  -q  •  . (q)  >  -q  . 

n*i  n,n-i 


(2.17) 


as  well  as 


~  -  >  0,  all  1 


(2. IB) 


Proof.  By  (2.1) 


K  (X jq)  »  0  Iff  H  (q"/X>q>  -  0  . 
n  n 


Since  we  have  ordered  y  (q)  as  In  (2.10),  (2.15)  follows.  (2.16)  Is  a  special  case  of 

n,i 


(2.15).  Fran  (2.15)  and  (2.12)  we  find 


«"  *  Vi(q,Vn-i(q>  51  <jn’2iVn-i(q)  ' 


which  implies  (2.17).  Finally,  combining  (2.15)  and  (2.7)  we  obtain  (2.18). 


Tl  "  1 

Theorem  2.2.  Let  q  >  1.  Then  for  i  *  1,2,...,  [ — r — J 


n-i  ,  ,  .  „  n-i 

"c1q  >'n,i(q>  *  _c2q 


(2.19) 


—  q*  <  y  , (q)  <  -  ~ —  q* 
Cj  M  n,n-i  M  c  H 


(2.20) 


the  constants  c,,  c2  do  not  depend  on  q  and  1,  and 


c  -  ly  ,(1))  , 

1  ll|  I 


(2.21) 


1  <  c2  < 


n  ♦  1 


n  -  1 
n  +  2 


,  n  odd 


(2.22) 


n  -  2 


,  n  even 


Proof.  It  is  enough  to  prove  (2.19),  since  (2.20)  follows  from  it  by  (2.15).  Observe  from 


n-1 


(2.11)  that  y  (q)  >  y  (1)q  .  Then  by  Lemma  2.1 


Vi<q)  >  ‘TT"  *  'in,t(1)q 


n-i 


_ciq 


n-i 


and  the  left  inequality  is  proved. 

Since  m  ,(«)/qn  *  is  a  continuous  function  on  [1,“t  and  satisfies  (2.17),  while 
n,l 

by  Theorem  2.4 


U 

lim  - 

q+»  g 


Uti 

n-i 


there  obviously  exists  a  constant  1  <  <  min  ~—£  —  independent  q  such  that  the  right 

v 

inequality  of  (2.19)  holds.  In  particular  note  that  1  <  c2  <  ^  for  n  «  2k  -  1.  ■ 

Theorem  2.2  bounds  p  , (q)  as  functions  of  q.  However,  it  is  of  interest  also  to 
n,l 

ask  the  opposite  question:  suppose  v  ,(q)  “  U  ,  .(g).  What  can  we  say  about  q,  q?  We 

n,i  n,i-i 

believe  its  answer  is  beautiful  enouqh  to  deserve  its  place  in  the  paper. 

Theorem  2.3.  There  exists  a  constant,  const  <  1  so  that,  for  any  q,  q  or  i, 

U  .  (q)  -  U  ,  ,(q)  >  (2.23) 

n,i  n#i-i 

implies 

<  const  <  1  . 

_  q 

Proof.  Let  q,  q  satisfy  (2.23)  for  some  i.  Then  (2.18)  gives  us 


n  n  n  u  (q) 

(2)  <  (?)  »  q  :*  a  solution  of  (f)  -  —-2i — 7^7  -?  P,(q)  . 


Vi-1,q)  ‘  1 


The  function  p^q)  is  a  continuous  function  of  q,  and  by  Lemma  2.1  and  (2.9) 


q  >  1  :  P^(q)  <  ~  , 

.  .  Vn-l+1,V<’) 

q  < 1 5  pi(e')  •  -? — "nTqT 

n  ,n— 1 


<  q 


Thus  Pi(0'*-)  »  Pi(»)  “  0.  Clearly  we  find 


7- 


const  -  max  max  p. (q)  <  1  . 
i  <J 


The  last  part  of  this  section  we  devote  to  the  asymptotic  behavior  of  u  (q)  as 

n#i 


Theorem  2.4.  For  1*  1,2,...,  f 


n_=JLl 

2  lf 


n  -  1  n-i  n-1-1  .  n-2-i. 

<<<?>  - t —  <T  -  C  .q  +  0(q  ) 

n#i  l  n,i 


(2.24) 


1  1  /  1  ^  1-1  1-2. 
•‘n.n-i^’  “  -  «  +  cn,i<?  +  0(ej  >  * 


(2.25) 


Here 


0  <  c 


n,i  i2(l  +  1 ) (n  -  1)(n  -1+1) 


t(n  -  21)4  +  (61  -  1 ) (n  -  21 ) 3 


+  41(31  -  1 ) (n  -  2i)2  +  4i2 ( 21  -  1)(n  -  21)]  . 


(2.26) 


In  particular. 


VKZ  1) 


2k-1,k-1  fc(k  _  1}2(k  + 


(2.27) 


Proof,  By  (2.15)  it  is  enough  to  prove  (2.24).  Let  X  -  y  (q).  By  (2,8)  and  (2.15) 
^  n#i 


lim  -  0,  all  1 

q+«  qn 


Thus  for  some  a  *  0  and  some  r  >  0 

A  -  <>q"'r  +  0qt''r_1  +  0(qn”r”2 )  .  (2.28) 

Since  the  coefficients  ain ^ <?)  are  polynomials  in  q,  and  after  a  proper  normalization 
in  1/q,  r  is  an  integer.  From  Corollary  2.1  we  conclude  that  as  q  ♦  * 


n-1 

I 

1-0 


n,l 


,  o1  .  V  (n+i)  (n-l-i)/2+(n-r)i.  (0)  1  -1 . 

(q)A  *  l  q  (a  a  +q  (i 

i-0  n,i 


d)  i  (0)  i- 
, a  +a  . ia 
n,i  n,i 


6)) 


(2.29) 


An  inspection  of  the  exponent 


8 


shows 


td, r)  (n  +  il(n  -  t  -  l)/2  +  (n  -  r)i 


»(2<n  -*■>  -  1  ♦  i.r)  -  ,  i<r)  , 


A  *(i.r)  s-  i(i(i  +  I.r)  -  Ki(i,r) 


r  -  1  -  1  . 


Since  6^i|i(n  -  r  -  I.r)  «  0,  the  leading  power  of  a  occurs  in  the  terms  i  =  n 


1 , 


n  -  r.  Thus  (2.29)  can  vanish  precisely  for  r  «  1,2, 


,n  -  1  as  q 


and  we 


conclude  from  (2.10): 
complete  the  proof. 


n-i. 


M  (q)  *  0(q  ).  By  using  (2.6)  it  is  now  straightforward  to 

n,i 


i 


i 


-9- 


3.  Polynomial  Splines  on  a  Biinfinite  Geometric  Mesh 


To  start  more  generally,  let  t  :«  (t^)^  be  a  strictly  increasing  biinfinite 

sequence  with  t._  :=  lim  t.  ,  Is-  It  ,t  [.  Let  further 
i-»±“ 


mS  (I ) 

n  j  t 


{f  e 


cn'2(i)  n 


L_(I) 


is  a  polynomial  of  degree  <  n} 


be  the  normed  linear  space  of  polynomial  splines  of  order  n  with  the  breakpoint 

sequence  and  the  norm  If  I  :*  sup  |f(x)|.  Let  r,k  e  N  be  given  integers, 

xel 

0  <  r  <  2k,  0  <  k.  Consider  the  map 

\  !  mS2k-r,t(I)  +  !  f  ~  (*i.rf)- 

associated  with  interpolation  conditions 


*1,0 


f(t.  ), 


*i,r 


/  M  f,  r  >  0  . 
I  i,r 


(3.1) 


Here/  as  usual  the  B-splines  of  order  k  with  knots  t^  are  defined  by 

k-1 

Mik,xl  k(ti'ti+r”*,ti+k] (*  "  x>+  * 

Nik  !*  i  {ti+k  -  VMik  * 

The  interpolation  problem:  for  given  b  :»  (bj *  im,  find  f  e  s2)c_r>t*I>  such  that 

Rrf  -  b 

is  by  [21  correct,  if  Rr  is  invertible,  i.e.  the  Gramian  (totally  positive)  matrix 

Gr  (*i,rNj,2k-r)i,j— » 

is  boundedly  invertible. 

i  +■> 

Let  us  restrict  ourselves  now  to  a  particular  geometric  knot  sequence  t  :*  (q 
for  some  q  t  1 0 ,°»t .  In  this  case  the  matrix  is  a  Toeplitz  matrix  and  is  boundedly 
invertible  iff  the  characteristic  polynomial 


5  >  Yr’W 


(3.2) 


has  no  zero  on  the  unit  circle  |1|  -  1,  or  since  Gr  is  totally  positive,  at  X  -  -1. 
The  case  r  »  0  is  treated  in  [6] ,  where  it  is  proved  that 


-10- 


WX!<> 


2k -2 

I 

3-o 


I  *3* 


2k-1,0Nj,2k 


X“  ‘  '  \  XVo»},2k 


and  from  properties  of  the  generalized  Euler-Frobenius  polynomial  determined  when  R0  is 
invertible.  A  nice  argument  shown  to  us  by  de  Boor  [3]  leads  to  the  conclusion: 

The  characteristic  polynomial  (3.2)  has  -1  as  a  zero  iff 

n2k-1<-q’t‘r,<?)  “  0 

for  any  r,  0  <  r  <  2k  -  1.  A  recent  result  of  Hollig  [5)  states 


.-1 


■Gr  -  hr(q) 


Wq  'q> 


n2k.1(-q  «q> 


He  proves  that  hr(q)  is  bounded  independently  of  q  and  Gr,  r  *  k  -  1  ,k  is  not 
boundedly  invertible  for  at  least  one  q  e  f1,»[.  We  give  here  an  alternative  proof 
simply  rereading  Theorem  2.2.  By  Theorem  2.1  we  can  restrict  to  the  case  0  <  r  <  k 
The  equation 

Vq)  !*  ni.k.r<q)  !“  w2k-1,i(q)/qr  ■  -1 

has  (at  least  one)  solution  q  e  ]0,»(  exactly  for  r+l<i<2k-2-r.  Put 

cr  »”  (qlq  is  a  solution  of  (3.5)} 

and  |Qr1  i*  number  of  elements  in  Qr.  Choose  r+1<i<k-1.  Then 

2k-1-r-i  .  2k-1-r-i  .  , 

“C^q  <  (o)  <  -  c^q  ,  q  >  1  , 

~t+i  .  _  ,  -1  -r+1  ,  . 

-c,  q  <  n.(q)  <  -c  q  ,  q  <  1  • 


by 
-  1. 

(3.5) 


(3.6) 


If  q  >  1  obviously  there  is  no  solution  to  (3.5)  since  this  would  imply  i  >  2k  -  r.  In 
the  case  q  <  1  there  is  q  *  Cr  exactly  for  i  >  r  +  1.  Since 

ni(q)  -  -1  iff  n2k_1_1(l/q)  =  -1 

our  claim  is  confirmed. 

We  note  that  the  case  of  a  finite-partition  [1]  suggests  that  Rr,  r  *  k  -  1*V  is  not 
invertible  for  all  q,  but  as  already  pointed  out  in  (2]  the  same  proof  can  not  he  applied 
since  the  quotients 
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min 

i+1-r<j<i 


< 


r 


are  bounded  independently  of  q. 

Let  now  q  >  1.  From  (2.15)  we  get 


Vq) " 


^  -  U2k-1,l(q) 
i-1  *  *  ^(q) 


V-1  w  {q)  +  1 
|  |  — i - 

i-1  Wi<ql  "  1 


(3.7) 


with 

Vq)  *"  -<U2k-1,l<q)  +  ‘,2k-1,2k-1-i(q,)/(qC  +  qk'1)  •  (3‘8) 

From  Theorem  2.2  we  conclude 

—  k  v-i -i  -i  -v+i+i 

Wjfq)  *  <c2<»  +  «s2  q  )/2  w^q) 

and 


te-1  w  (q)  +  1 

h.  .(q)  I  I  ZT -  <  h.  .(q)  <  h.  .  (o) 

*  1-1  w  (q)  -  1  * 


k-1  w. (o)  +  1 

I  |  1 

i-1  2i(q>  "  1  ’ 


Sine.  is  decreasing  as  a  function  of  q,  this  suggests  that  h)c_1(q)  is  too. 

However,  we  succeeded  in  proving  this  only  as  q  ♦  as  a  consequence  of  Theorem  2.4  and 
(3.7),  (3.8). 

Theorem  3.5.  For  0<r<2k-1,  the  Gramian  matrix  is  not  boundedly  invertible  for 

q  *  and  ICrl  “  IC2k-1-r*  *  2(k  In  particular, 

lBkl  -  end  th*  no*™  for  <1  *  f1,*f  satisfies 


also  as  q  +  - 


h^tq)  <  h^tq)  <  h^q)  < 


fC2  *  1} 

c,  -  1 1 


2<k-1) 


hy(q)  -  (2k  -  1)(1  +  4(k-2)/(k+1)q-1  +  0(q’2))  . 


(3.9) 


(3.10) 
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